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ABSTRACT 



FUSE spectra were obtained of the supersoft X-ray binary RX J0019. 8+2156 
(QR And) during 16 consecutive spacecraft orbits, covering the binary orbit 
(P=15.85 hr) with about 0.2 phase overlap. The spectrum is dominated by 
strong H2 absorption (column density ^10^*^ g cm~^), which appears at the 
velocity different from other interstellar absorption lines and may be partially 
circumbinary. This absorption makes study of spectral features from the binary 
system difficult. The only well-detected emission lines are He II, lOSSA, and 
O VI, IO32A, (the other line of the O VI doublet, 1037A, is largely obscured by 
strong H2 absorption). The O VI shows a P Cygni profile that varies in velocity 
and strength with binary phase. We compare this with similar changes seen in 
Balmer line profiles. We extract the FUV light curve and compare it with the 
optical light curve. There is an eclipse in both wavelength regions, but the FUV 
minimum lasts much longer, well beyond the visible light egress. The FUV 
results are discussed in connection with the binary model and mass fiows within 
the system. 

Subject headings: ultraviolet: stars - (stars:) binaries: close - X-rays: binaries - 
stars: individual (RX J0019.8+2156 = QR And) - ISM: jets and outflows 
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1. Introduction 

RX J0019. 8+2156 (hereafter called by its variable star designation, QR And) is one of 
the few identified Galactic supersoft X-ray binaries, although numerous examples have been 
found in the Magellanic Clouds and M31 (see Greiner 1996). The system was first noted in 
the ROSAT All-Sky Survey and described by Beuermann et al. (1995). They found this 
12th magnitude system to show ~0.5 mag eclipses and radial velocity variations which fit 
an orbital period of 15.85 hrs. Archival photographic plates show the orbital period has 
been stable for over 40 years (Greiner & Wenzel 1995), although there are irregular changes 
of a few tenths of a magnitude in the overall brightness level of the system (e.g. Cowley et 
al. 1998, Deufel et al. 1999) 

The system is at Galactic latitude —40°, and hence shows little evidence of interstellar 
reddening, with mean U — B ^ —0.7 (Beuermann et al. 1995) and quite weak interstellar 
lines of Ca II and Na I (McGrath et al. 2001). This, with the mean absolute magnitude of 
supersoft binary systems, places its distance at ~2 kpc, so it lies well out of the Galactic 
plane. 

Optical spectra show very strong emission lines of He II (virtually all of the He II lines 
in the Multiplet Table are identified). Balmer absorption lines are present through about 
half of the orbital cycle and indicate a strong outfiow of material with velocities of ~ —200 
to —600 km s~^. Like several other supersoft X-ray binaries, QR And also shows bi-polar 
emission-line jets (velocity ~ ±800 km s~^) in He II (4686A) and the lower Balmer lines 
(e.g. Becker et al. 1998, Cowley et al. 1998, Deufel et al. 1999). These jet features come 
and go on timescales of weeks to months. Since the Balmer outfiow is phase-dependent (and 
persistent?), this is likely a constant disk fiow that is not strongly related to the jet activity. 

There has been some disagreement about the appropriate model for QR And, with 
various authors deriving quite different values for the orbital inclination, and hence the 
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stellar masses. Beuermann et al. (1995) concluded that the inchnation must be quite low 
{i ~ 20°), based partially on the assumptions that the masses be similar and consistent 
with a white dwarf. On the other hand, Meyer-Hofmeister et al. (1998) analyzed the light 
curve and presented models which indicate an inclination near i = 55°, while Tomov et al. 
(1998) suggested an even larger value {i ~ 79°) also based on the light curve. The adopted 
inclination is important since it ultimately confines the range of possible masses, and hence 
the model and evolution of the system. /^From an analysis of the emission-line velocities 
and those of the bi-polar jets, Becker et al. (1998) found that the jet lines show the same 
velocity amplitude and phasing as He II, 4686A, and hence concluded that the He II lines 
must reveal the motion of the compact star where the jets originate. Further, from analysis 
of the outfiow line profiles they found that the orbital inclination must lie in the range 
35° < i < 60°. For the purpose of our discussion below, we have adopted these conclusions 
of Becker et al. and assume a model similar to that discussed by Meyer-Hofmeister et al. 

O VI emission (3811, 3835, 5290 A) is found in all supersoft X-ray binaries and is 
potentially interesting as its high ionization may indicate an origin close to the accretion 
disk center. However, study of these lines is compromised by blends with He II at 3813A 
and 3834A, and with [Fe XIV] at 5303A. This suggested that it would be particularly 
interesting to observe the supersoft systems in the far ultraviolet with FUSE since the O VI 
resonance lines (1032, 1037A) lie within its wavelength range. 

2. FUSE Observations 

An overall description of the FUSE data is given by Sahnow et al. (2000). Our FUSE 
observations were carried out during 16 successive spacecraft orbits on 2000 July 27-28. 
The large apertures (30" x 30") were used. No unusual anomalies were apparent during 
the observation sequence. The observations were taken in time-tag mode, and the 16 
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separate spectra cover slightly more than one binary orbit. Table 1 lists the details of the 
observations. 

The spectra were processed by the CALFUSE pipeline and retrieved from the archive. 
Later, the data were reprocessed using an updated version (1.8.7) of the pipeline processes. 
The results were the same to within about 1%. 

Surprisingly, the spectra were found to be heavily absorbed by molecular hydrogen 
(H2). The amount of absorption is characteristic of more reddened targets (Eb-v =0.5 or 
more), whereas, as noted above, QR And has quite blue colors. The spectra also contain 
airglow emissions (see Feldman et al. 2001) which are much stronger on the daylight side of 
the FUSE orbit. Spectra 6, 8, 10, and 12 were taken almost entirely during orbital night, 
and thus they provide a way to assess the airglow contamination of the other spectra. 
Fortunately, the spectral region near the O VI 1032A emission line is free of airglow lines, 
so that our analysis of this region could be performed on data from both orbital night and 
day without editing. 

3. Spectral Energy Distribution and H2 Absorption 

In order to assess the contamination due to H2 absorption, we used simple absorption 
models with a range of column densities. We also derived an empirical absorption model by 
using a median spectrum of the region of interest around the O VI resonance lines. The 
1032A line is free of strong H2 absorptions, but the 1037A line is almost completely overlaid 
by saturated absorption. 

Figure 1 shows the average spectrum with different amounts of smoothing and airglow 
removal, as well as the H2 model from our grid that fits it most closely. The H2 absorption 
strength in the FUSE wavelength range is loosely correlated with Eb-v and was compared 
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with several BO stars with E^_y ranging from 0.15 to 0.5 magnitudes which had been 
observed with FUSE. The values for QR And correspond to values of Eb-v ~0.5. However, 
Beuermann et al. (1995) estimated a value of E^^y ~ 0.10 for QR And, based on the 
strength of its 2200A absorption feature. The colors of QR And {U — B = —0.7 and 
B — V = 0.01) also show the system to be very blue, so the E^.y cannot be 0.5 and the 
H2 absorption strength is unusually high for the reddening. Since the system is losing 
matter via disk and jet outflow, it is worth considering whether some of the H2 absorption 
is due to material surrounding the binary. We estimate the column density causing the 
H2 absorption to be ^10^*^ cm^^. Beuermann et al. (1995) estimated an absorbing column 
of 4 X 10^° H- atoms cm~^ from ROSAT data for a combination of Galactic foreground 
absorption plus any possible intrinsic absorption. 

To determine if the H2 absorptions are interstellar or circumstellar, we compared their 
velocities with Ca II and Na I interstellar lines and with QR And's systemic velocity, both 
derived from optical spectra. From the FUSE spectra we measure the heliocentric velocity 
of H2 features near the O VI lines to be —16 ± 2.5 km s~^, with no variation from spectrum 
to spectrum. The absolute wavelength precision depends on placing the star accurately in 
the (guided) LiFl channel aperture. A precision of 1 arcsec is expected, which corresponds 
to ±3.4 km s~^. To further check the FUSE wavelength zero-point (possibly due to grating 
movement uncertainty) we measured the Ly/3 airglow wavelength. This tracks the orbital 
velocity of FUSE along the boresight to within 2 km s^^, so we conclude that velocities are 
correct as measured. 

An MMT spectrum obtained 1995 Oct. 12 shows the mean velocity from interstellar 
absorption lines of Ca II and Na I to be +4 ± 4 km s^^. McGrath et al. (2001) found the 
systemic velocity of QR And to be ~ —41 ±3 km s~^ (Beuermann et al. (1995) gave —59 ±2 
km s~^ from lower resolution data). The measured H2 velocity of —16 km s~^ lies between 
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the ISM and binary velocities. In view of the unusual H2 strength, it thus seems likely that 
it is partly associated with the binary system - perhaps formed during a PN episode and 
well-removed spatially from the current ionizing radiation from supersoft X-ray activity. 

4. Far Ultraviolet and Optical Light Curves 

The FUSE data were used to construct the FUV light curves shown in Figure 2. FUSE 
consists of four separate telescopes that are kept co-aligned by compensating for modelled 
mechanical distortions around the orbit (e.g. Moos et al. 2000). Only one telescope (LiF-1) 
is used for tracking. The extraction of signal is thus uncertain in the other three telescopes 
if the signal is weak and the position of the target in the (30") aperture drifts by a great 
deal. The LiFlA (longer wavelengths) are also subject to a detector wire shadow (the 
'worm') that may affect photometric accuracy. However, the LiF2A (longer wavelength: 
1090-1 I80A) channel has little drift and the highest signal level. Thus, we used the fluxes 
from the LiFlA (990-1080A) channel and the LiF2A channel to derive two independent 
light curves. The continuum fluxes were estimated by summing these spectral regions after 
removal of airglow lines and interpolating across the deep H2 absorptions. These are shown 
as a function of phase in Fig. 2, displayed as magnitudes. There is good correspondence 
between the two ultraviolet regions. The count rates for all channels are steady and show 
no variations that might indicate the star was drifting out of the aperture. The spectral 
images show that the extraction windows appear to cover the data well. In the binary 
phases where there are overlapping FUSE observations (phases 0.5-0.7), the fluxes are in 
good agreement. 

The FUSE fluxes were rechecked by one of us (AWF) on the raw data using the current 
version of CALFUSE (1.8.7). The calibrated spectra include the total counts associated 
with each pixel (i.e., the sum across the width of the spectrum in the spatial direction). A 



light curve from the total signal over the interval 1 102-1 133A was derived from the LiFlB 
channel. This region does not have any strong airglow lines and avoids the LiFlB channel 
'worm'. The result is essentially the same as the 1090-1180A LiF2A curve shown in Fig. 
2. To check CALFUSE further, spectra for 5-minute intervals were extracted from the raw 
data, and count rates were determined over the same 1102-1 133 A wavelength region. These 
light curves are also very similar to the FUV light curves shown in Fig. 2. 

Comparing the two FUV light curves in Fig. 2, we find that overall amplitude is larger 
by ~ 0.03 magnitudes in the longer wavelength channel. This is likely to be caused by 
H2 line blanketing that lowered the continuum in the shorter wavelength region where the 
absorptions are more crowded (see Fig. 1), rather than being a real color change. We note 
that the visible light shows almost no color changes through eclipse (e.g. Cowley et al. 
1998, Deufel et al. 1999). 

For comparison of the optical and far-ultraviolet fiuxes. Fig. 2 shows two l^-band light 
curves together with the two derived from the FUSE data. The nearly complete optical 
light curve (small filled circles) was obtained in 1995 September and is fully described by 
McGrath et al. (2001). It is plotted with 9-point smoothing applied. In the unbinned 
optical data, fiickering is present throughout the orbital cycle, and quasi-periodic variations 
(P ~1.8 hours) are prominent at orbital phases 0.2 to 0.6. The second optical light curve 
(small open circles) shows V data taken on 2000 July 28 by one of us (PCS) at Braeside 
Observatory in Flagstaff, Arizona, with monitoring beginning just hours after completion 
of the FUSE observations. The orbital phase coverage is rather limited, but it includes 
most of primary eclipse. The short-term fluctuations are of much smaller amplitude than in 
the 1995 data. In 2000, QR And was ~0.1 magnitude brighter than in the 1995 data, but 
variations in the mean brightness level are well-known in this system (cf. Will & Barwig 
1996, Matsumoto 1996, Cowley et al. 1998, Deufel et al. 1999). The optical data shown 
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here illustrate the typical behavior of this source. 

There are very obvious differences between the optical and FUV light curves in the 
phases following minimum light. While the ingress is similar in both wavelength regions, 
the FUV light curve shows a minimum that lasts from phase 0.95 through 0.25, much 
longer than the minimum observed in optical light. Thus, we appear to have an eclipse 
of the hottest regions of the disk that extends well beyond that of the cooler regions seen 
in optical wavelengths. If this extended minimum is real, it could indicate important disk 
structures in the system which presumably effect the FUV radiation. However, we would 
really like to observe another minimum to be certain there was no FUSE instrumental 
problem. The 'secondary eclipse' dip present in the visible light curve is not seen in the 
FUV, but a small dip does occur at phase ~0.6. 

If the FUV extended eclipse represents a permanent feature, we might expect to find 
a gradual change in the light curve as one observes at progressively shorter wavelengths. 
The best we can do at the present time is to examine other light curves available in the 
literature. Beuermann et al. (1995) show light curves at 1400-1500A and X-ray data 
from RO S AT-PSPC . Replotting these data with a more extended vertical scale we find 
both the duration of minimum and the phases of ingress/egress to be consistent with our 
optical light curve and that of Will & Earwig (1996). Similarly the lUE data shown by 
Gaensicke, Beuermann, & de Martino (1996) for 3000A and 1250A when replotted also 
show similar phases of ingress/egress and duration of minimum flux. We were unable to 
find any indication from published ultraviolet data of other extended-duration eclipses 
such as seen in the FUSE data. However, Gaensicke et al. (2000) show ROSAT-B.KI 
X-ray data obtained between 1997 Dec. 2 and 1998 Jan. 10 in which there is no evidence 
of any orbital variation. Thus, changes do occur in QR And which severely modify the 
light curve. Furthermore, Meyer-Hofmeister et al. (1998) found that important changes in 
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the disk-rim structures can occur on timescales as short as the orbital period. However, 
without additional FUSE data, we don't know if the extended minimum we observed in the 
FUV light curve was a transient event or due to a more permanent structure in the binary 
system. 

5. FUV Spectrum 

The spectra were examined carefully for the presence of high ionization lines in the 
FUSE range. While the H2 contamination makes this difficult, we find no evidence of 
N III, C III, S IV, Si V, or P V. However, O VI is clearly present and has a complex and 
variable line profile that we discuss in detail below. The only other line we detect is He II 
IO85A. This line is close to two airglow lines of N II, but both the median of all spectra 
and the sum of the night spectra show there is a broad He II emission line present in 
addition to the narrow and variable airglow lines. Due to the airglow contamination and 
the weakness of He II, we are unable to study its behavior as a function of binary orbital 
phase. Nevertheless, it is clear that the line does not have the P Cygni absorption profile 
found in the O VI lines. This is consistent with the purely emission profiles of the visible 
Pickering He II fines reported by McGrath et al. (2001). We note that in ORFEUS II 
observations (900-1200A) no stellar emission lines were detected (Gaensicke et al. 2000). 

In the optical region, lines of C HI and C IV (and possibly N HI and N V) are weakly 
present (McGrath et al. 2001), but lUE spectra of QR And show no evidence of the 
resonance lines of C IV at 1550A (Beuermann et al. 1995). The only emission feature in 
the lUE data is He II 1640 A. In the FUSE spectra, C HI at 11 75 A is not detected, even 
though it is clear of H2 and airglow. We set an upper limit on its equivalent width of less 
than O.lA. It appears that carbon, and possibly nitrogen, have low abundances in QR And. 
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5.1. O VI 



The strongest emission line in the FUSE range is O VI, 1032A. This whole feature is 
fortunately free of airglow or major H2 absorption (there are some weak narrow absorptions 
but they are comparable with the noise levels in individual spectra and very much narrower 
and weaker than the stellar lines). The line has a P Cygni profile which changes around the 
binary orbit, as shown in Figure 3. By subtracting an H2 model and looking at differences 
from the median of all spectra, we can also see the edges of the absorption and emission of 
the O VI IO37A line, which gives us added confidence that the changes we see in the 1032A 
line are real. 

The extensive spectral coverage of the optical region by McGrath et al. (2001) shows 
that the lines near Balmer wavelengths have P Cygni profiles that change markedly around 
the orbit. (Note that the emission portion of these lines is a blend of H+He II, as one can 
tell by comparing their strengths with the non-blended, adjacent He II Pickering emission 
lines. However, the variable-strength absorption component is due only to H.) Fig. 3 shows 
the changes in both the FUSE O VI line and the optical 'Balmer' lines during similar 
phases (computed using the optical ephemeris of McGrath et al.) The FUSE data are 
smoothed by 91 points, which makes the resolution comparable with the optical, without 
losing the significant changes in the broad line profiles. In order to compare the O VI and 
'Balmer' lines, the optical spectra were averaged from numerous spectra within the phase 
bins specified, but were not smoothed. The profiles in Fig. 3 were normalized to the local 
continuum, and measurements were made of the centroids and fluxes of the emission and 
absorption components as shown in Figure 4. The values plotted in Figure 4 have measured 
total scatter of ±20% for the EW values, and velocities are reproducible to ~50 km s~^, 
taking into account uncertainties in continuum fitting. 

The absorption outflow velocities are similar for all three lines shown in Fig. 4, but the 
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variations in absorption strength differ. O VI has absorption at all phases, while little or no 
Balmer absorption is present from phases 0.7 to 0.0. The O VI and H absorption strengths 
thus show no binary phase correlation in Fig. 4, while the O VI emission and absorption 
do follow the phase variations of the Balmer lines, albeit with a smaller amplitude. 
This suggests that there may be significant azimuthal variations of the ionization of the 
outflowing material. In this connection we note that FUSE data show anomalous O VI 
presence in the outer regions of OB star stellar winds (Crowther, Bianchi: in preparation) 
whose ionization may be connected with X-ray flux. Further work is needed to understand 
these phenomena properly. 

If the velocity of the O VI emission is due to orbital motion of the compact star, 
maximum velocity should occur at photometric phase 0.25. However, the maximum occurs 
~0.1P earlier, almost certainly due to the very strong P Cyg absorption present between 
phases 0.0 and 0.2. The O VI velocities are very similar to those found by McGrath et al. 
for the blended H+He II emission lines. These authors suggest that only the pure He II 
lines (i.e. those with no trace of overlying H absorption) reveal the orbital motion of the 
compact star, the same conclusion as was reached by Becker et al. (1998). The variable 
strength of the hydrogen absorption modifies the emission-line profiles, hence distorting 
their velocity curves. 

Changes in the line profiles are shown in Figure 5. Differences from the phase bin 
showing emission-only 'Balmer' lines are plotted, so the figure shows the changes in the 
absorption with phase. The O VI absorption is more complex and is seen at all binary 
phases. It is weakest during the low parts of the FUV light curve, suggesting that the O VI 
absorption is associated with the FUV continuum source. 
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6. Discussion and Summary 

The optical spectra of McGrath et al. show that there are significant outfiows from the 
system, with strong azimuthal variations. Generally speaking, the outfiows have velocities 
of several hundred km s^^ and are seen most strongly through mid-eclipse and egress, but 
not seen during eclipse ingress. In the FUSE data, the O VI resonance lines show an outfiow 
of very hot gas at all phases, but with velocities similar, although somewhat lower, to the 
Balmer absorption. 

The emission and absorption features in O VI and H+He II show similar velocity 
changes with binary phase. These are not the same as the pure He II emission velocity 
curves which appear to trace the binary orbital motion (e.g. Becker et al. 1998; McGrath 
et al. 2001). Thus, the wind from the binary has azimuthal velocity variations. The 
relative absorption strength also changes with binary phase, which may indicate azimuthal 
ionization changes. The production of O VI in OB stars is anomalous and may have some 
connection with X-ray ionization which could also be occurring in QR And. 

We find that the high H2 absorption column, which is anomalous for the unreddened 
system color, has a different velocity from the atomic interstellar lines, lying between the 
ISM and assumed binary systemic velocity. We suggest that it may partially arise in a 
circumbinary location, possibly formed during the presumed PN episode that occurred 
prior to the present evolutionary stage. We note that circumbinary II2 is detected around 
some planetary nebulae and that other degenerate binaries we are studying with FUSE do 
not show appreciable H2 absorption. 

Finally, the very odd FUV light curve shows an extended eclipse minimum lasting 
about one third of the orbit, although only one eclipse was observed by FUSE. While it may 
be tempting to connect the extended FUV eclipse with variable II2 absorption, the FUV 
light curve was derived from the longer FUSE wavelengths which are free of H2 absorption. 
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Furthermore, we have determined that there is no change in the column density of H2 or 
its velocity during these phases. The difference between the visible and FUV eclipse might 
be explained by a UV-bright spot on the rim of the disk, where a mass-transfer stream 
impacts it (e.g. models by Meyer- Hofmeister et al. 1998). The off-axis location of such a 
spot could give rise to an extended eclipse in the FUV. Such a model would imply that the 
effect becomes more visible as the wavelength of observation decreases from visible to FUV. 
Our examination of archival photometry at UV and X-ray wavelengths suggests that the 
light curve may show large long-term changes. 

If the O VI line profile (both emission and absorption) arises in part (or entirely) at 
this UV-bright location, the observed phasing of the O VI radial velocities is qualitatively 
as expected. One can also imagine different outflow velocities for O VI, as observed, if 
the line arises primarily from this location rather than from the entire disk. The lack of 
hydrogen absorption from phases 0.65 to 0.95 might occur if the region around this hot spot 
is highly ionized. 

APC and PCS gratefully acknowledge support through NASA grant NAG5-8805. 
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Table 1. FUSE Observations of RX J0019.8+2156 



Orbit Photometric phase" 



HJD (mid) 


Exposure 


2451700+ 


(sec) 


52.922 


2481 


52.993 


2191 


53.066 


1871 


53.135 


1590 


53.206 


1310 


53.246 


418 


53.277 


1160 


53.318 


850 


53.347 


920 


53.389 


1220 


53.418 


701 


53.461 


1721 


53.4784 


1320 


53.541 


3461 


53.610 


3461 


53.679 


3301 



1 0.475 

2 0.584 

3 0.694 

4 0.799 

5 0.907 

6 0.966 

7 0.013 

8 0.075 

9 0.120 

10 0.184 

11 0.227 

12 0.293 

13 0.327 

14 0.414 

15 0.519 

16 0.622 



^ephemeris: To=HJD 2449987.8459(16) + Ex0.6604645(14) days (McGrath et al. 2001) 
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Fig. 1. — (upper two panels) Average FUSE spectrum of QR And with different amounts 
of smoothing and airglow removal, (bottom two panels) For comparison, the H2 model 
absorption is shown with two different smoothings. 

Fig. 2. — Light curves for QR And. The small open circles (upper curve) are V data 
obtained on UT 2000 July 28, only a few hours after the FUSE observations ended. The 
small filled circles are the V data of McGrath et al. (2001). It is known that this binary 
shows variations of a few tenths of a magnitude in its mean brightness level (e.g. Cowley 
et al. 1998, Deufel et al. 1999). The lower two curves are the FUSE light curves extracted 
from different wavelength regions, as indicated. Phases are computed using the ephemeris 
given by McGrath et al. 

Fig. 3. — Line profiles of O VI, 1032A, binned through the orbital phases indicated. These 
are compared with line profiles of H/3 and H7 in the same phase bins (from the data of 
McGrath et al.). The FUSE data have been heavily smoothed to spectral resolution ~ lA; 
the optical data had an original resolution of ~ SA, but many spectra have been averaged 
together. The vertical scales are in units of the local continuum. 

Fig. 4. — Velocity and equivalent width measurements from FUSE and optical spectra 
binned as in Fig. 3. The measurements are made with respect to the continuum levels 
plotted in Fig. 3. In the upper panel, the mean of H7 and H/? is plotted to simplify the 
diagram. The Balmer velocities are not corrected for the He II blending, which would be 
—60 km s~^ for equal contributions of H and He II to the emission. 

Fig. 5. — Similar to Fig. 3, but showing changes in profiles from the 0.67-0.95 phase bin, 
when Balmer absorption is not seen. The phase 0.67-0.95 profiles are shown in the lowest 
plots. This figure illustrates the changing absorption components with binary phase. 
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